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Abstract 

Motivated by the Dirac idea that fundamental constant are dynam- 
ical variables and by conjectures on quantum structure of spacetime at 
small distances, we consider the possibility that Planck constant h is 
a time depending quantity, undergoing random gaussian fluctuations 
around its measured constant mean value, with variance and a 
typical corrélation timescale At. We consider the case of propagation 
of a free particle and a one-dimensional harmonie oscillator cohér¬ 
ent State, and show that the time évolution in both cases is different 
from the standard behaviour. Finally, we discuss how interferometric 
experiments or exploiting cohérent electromagnetic fields in a cavity 
may put effective bounds on the value of r = 
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1 Introduction 


The idea that fundamental constants in physics may be spacetinie varying 
parameters was first suggested by Dirac with his “Large Numbers Hypoth- 
esis” [^. Since then, many efforts hâve been devoted to embed this idea in 
cohérent theoretical frameworks. One possibility is that these constants may 
be effective quantities. This happens, for example, in higher-dimensional thé¬ 
ories, after réduction of coupling constants to the four-dimensional subspace. 
For a critical summary of this subject see [^[^. 

The most studied examples of varying fundamental constant are perhaps 
the Newton gravitational constant Gn and the hue structure coupling a. In 
the latter case there is quite a large number of experimental constraints on 
its time évolution, spanning several order of magnitude in redshift, or cosmo- 
logical time, from laboratory experiments, geophysical tests as the Oklo phe- 
nomenon, the anisotropies of Cosmic Microwave Background, till the early 
epoch of Big Bang Nucleosynthesis, see and references therein. Spatial 
and time variations of a hâve been also constrained by observations of ab¬ 
sorption spectral lines in intervening clouds along the line of sight of distant 
Quasars at redshift s ^ ~ 2 — 3, |4]-[^. 

Evidences for dynamical features of fundamental constants would hâve 
of course, an enormous impact on the general picture of physics laws. The 
typical approach to this theoretical and experimental issue, as in the case of 
a just mentioned, is to look for signais of variation over long time period, 
billion years in the case one uses cosmological observables. The underlying 
idea is that the dynamical field which describes a particular fundamental 
constant is changing as a classical background held under the action of some 
effective potential which drives the held towards its minima. However, if this 
putative held T reaches a stable conhguration in very short time intervals, 
so that there is no further évolution in the long time régime, we could only 
see an ehect by strongly perturbing this conhguration and moving T away 
from the minimum. This can be achieved by exploiting the T coupling with 
matter helds and exciting the held by pumping energy into it, i.e. the typical 
accelerator approach to unveil new degrees of freedom. However, even if 
T seats at the classical potential minimum, we expect its value to undergo 
small random huctuations. In this case the question is how these huctuations 
might produce observable ehects in experiments. 

We recall that many theoretical arguments point ont that at short time 
and distance scales both quantum held theory and gravitational interactions 
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described by general relativity ought to be embedded in some form of quan¬ 
tum gravity, a theory not yet known. Whatever is this theory, it will deeply 
change the geometry of spacetime and phase space. Examples are théories 
in which spacetime is composed by branes, feature a noncommutative geom¬ 
etry or are characterized by the spin foam of loop quantum gravity, see, e.g. 
[ 9 , 10 and 11 and references therein. The rôle of a spacetime fuzziness and 


stochastic Lorentz invariance violation lias been recently studied in 
Since its introduction at the dawn of last century 
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13 , h lias been consid- 


ered one of the fundamental constant ruling physical phenomena. Its value 
lias been measured by different methods spanning the spectrum of black body 
radiation, kinetic energy of photoelectrons in photoelectric effect, and more 
recently, X-ray density crystal measurement, and watt balance. The latter 
two hâve recently achieved a relative précision of 1.2 parts in 10® 
actual accepted value is h = 6.626070040(81) ■ 10“^^ J 


14 . The 


15 . The précisé 


measurement of h plays a crucial rôle in establishing both a new SI and a 
revised fundamental physical constant System. 

The (reduced) Planck constant h enters into the dehnition of the Heisen- 
berg uncertainty principle, one of the pillars of the probabilistic nature of 
quantum mechanics. In fact, deformations of the quantum phase space 
hâve been put forward in the form of generalized uncertainty principle, see 
described by the commutation relation 


e.g. 


16 19 


[x*,pj] =ih {ô] + /j(x,p)) , 


( 1 ) 


where /* dépend upon sonie “small parameter”. This is particular case of a 
most general phase space commutation relation among position and momenta 
which also includes 


[x\x^'] = 

[p*,Pj] = iCij (2) 

where f,0C are such that Jacobi identities are satished. In general, these 
relations break Lorentz invariance, and a way to restore it to consider /, 6 
and C to be random variables with microscopie variations, so that for large 
distances and time intervals the invariance is restored. The commutations 
relations ([^ and (|^ mean that the phase space has a non canonical sym- 
plectic structure, but with a Darboux transformation it is always possible to 
express (locally) the structure as a nontrivial commutator between position 
and momenta, ail other commuters being zéro. Since the transformation 
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mixes positions with momenta this calls into question the meaning of the 
position and momentum. We will consider that at any given instant the ob¬ 
servables can be seen as generated by two sets of operators, x‘ and pj, with 
commutation relations as in Q. 

If one considers that corrections to the commutators are small and dépend 
on the background and fluctuations at some high energy scale, the net resuit 
is that the r.h.s. of commutation relations can be effectively described by a 
random function acting as an effective dynamical ûeld. This is tantamount 
to consider h as a fluctuating parameter. In this paper we will consider that 
the space fluctuations are integrated ont, and will concentrate on the time 
variations. Energy scale dependence of h induced by nontrivial commutation 
relations has been proposed in 1^. In 21 the possibility that h is a held 


whose expectation value can vanish at high températures was considered as 
a possibility to quantize gravity. 

A comment about the issue of dimensionality. Planck constant is not 
adimensional, and daims about its time dependence would be always re- 
lated to some particular combinations of h and other dimensionful quanti¬ 
fies, depending on the adopted experimental method. Only for adimensional 
combinations of fundamental constants, such as a, it is meaningful to speak 
about spacetime variation unambiguously. The variability of fundamental 
constants is a délicate issue, see e.g. 22 . We will see next that the stochas- 


tic nature of h is encoded in a dimensionful time parameter r and observable 
effects for light propagation or a harmonie oscillator show up through the 
dimensionless combination car, with u the light or oscillator frequency. 

Our aim in this paper is not to develop a specihc theoretical framework 
describing this scénario. We will only assume that there is no classical back¬ 
ground évolution, but the held nature of the Planck constant shows up in 
random huctuation around the constant measured mean value, which we will 
continue to dénoté by h. 

2 General formalism 


Random huctuations of h are introduced via an adimensional gaussian stochas- 
tic variable e{t), so that the effective Planck constant reads h{l + e{t)), with 


e{t) = 0, 


e{t)e{f) = T5{t — t') 


(3) 

(4) 
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Overline dénotés the mean over e probability distribution. The parameters 
are as follow: r = a‘^At, a is the variance and équation Q states that fluctu¬ 
ations are uncorrelated for time différences larger than a typical corrélation 
time At, analogous of the strength of fluctuation force in brownian motion 
dynamics in Stokes-Einstein theory for diffusion coefficient. These relations 
are time translation invariants. 

Planck constant fixes the commutator of canonically conjugated variables, 
position and momentum for non relativistic particles, and bosonic (fermionic) 
fleld and conjugate momentum fleld commutator (anticommutator). Con- 
sider a particle in one dimension. If h fluctuâtes, we hâve 

[x,p] = + £(f)). (5) 


The commutator with a Hamiltonian H gives the time évolution of operators 


dA 1 dA 


( 6 ) 


The flrst term on the r.h.s. now explicitly dépends on time via s{t), and 
the time évolution of A, the average over the probability distribution of 
shows non-standard behaviour, as we will see. Had we deflned the Poisson 
structure by normalizing the commutator to h{l + £{t)), ail effects cancel 
and we would be back to standard quantum mechanics. More interesting is 
the case of équation ([^, i.e. the fundamental Poisson bracket fluctuâtes. We 
investigate this possibility. 

Our approach is different from stochastic quantization 


23 , where a Sys¬ 


tem is considered in thermal contact with some external réservoir, and its 
quantum properties are obtained in the equilibrium limit. This is achieved 
by introducing an auxiliary fictious time. Equilibrium is reached in the large 
Actions time limit. In our case, Planck constant dépends upon physical time, 
leading to novel features. 

Quantum measurements is itself a probabilistic concept, and it only makes 
sense to compute observable averages over the probability distribution given 
by a wavefunction. We dénoté this averaging by (A). The relation between 
standard quantum mechanical average and average over e{t) distribution is 
discussed later. 

Commutation relation of équation ([^ is represented by deflning the action 
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of position and momentum on wavefunétions 'ip{x) 


x'0(a;) = A{t) xtfj^x) = , (7) 

pî/,(a;) = -ihB{t) ^ij{x) = B{t)poij{x), (8) 

with A{t)B{t) = 1 + e{t), and xq and po the canonical pair of standard 
quantum mechanics. A change of représentation A{t) —)■ A'{t) is a time 
dépendent dilation induced by the unitary operator 


U{t) = exp 



A'jt) 

A{t) 


(xoPo + PoXo) 


(9) 


Here we adopt the choice 

A{t) = B{t) = y/f + s{t ), (10) 


and assume that in this représentation équations of motion are given by (@. 
The rationale is that effects of any high energy scale theory producing a 
fluctuating h are expected to disturb phase-space elementary volumes. It 
seems unnatural that coordinates and momenta should be treated differently, 
and rescaled by different factors. 

Although in principle, ail values of e{t) are possible, we expect e{t) prob- 
ability distribution to be narrow, so the number of events corresponding to 
large values is exceedingly small. Since e{t) is constant in the time interval 
At, in a time T its value is randomly extracted T/At times. For T~ 13.8 
Gyr, the âge of the universe, and At > 5.4 ■ 10“^^ s, the Planck time, for a 
gaussian distributed e, the event number for £ < —1, i.e. imaginary A(t) 


T 

Ât 


— 1 


de 


1-00 V 2 


.2 

20-2 


TTCT^ 


T 


erfc 


v/2 




( 11 ) 


is smaller than unity provided cr < 0.06, a value which seems already too large 
to be acceptable. It is unlikely that in the whole history of the universe these 
large fluctuations were ever produced, and we can safely assume A(t) G M. 


3 Experimental strategies: two examples 

If Planck constant is randomly inconstant, the issue is to hnd experimental 
approaches to constrain the r, or to hnd signature of its dynamics. We discuss 
two possibilities, interferometric experiments and cohérent electromagnetic 
modes in a cavity. 
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3.1 Pree particles and long baseline inter fer omet rie ex- 
periments 

Consider a free partiale with mass m. Schrodinger équation reads 

= ^(l + ^W)Po^- (12) 


Since [po, H] = 0, fundamental solutions are still plane waves eigenfunc- 
tions of Po with eigenvalue po- In the représentation with M the xq 

spectrum 








(13) 


Using these solutions we can construct wave packets. Position and momen- 
tum operators act as follows 

= a/1 + e{t)x%lj{x,t) , (14) 

. (15) 

dx 

We want to evaluate position mean value and uncertainty versus time. As 
mentioned already, averaging should be performed both in the usual quan¬ 
tum sense, and over the e(t) probability distribution. We dénoté this double 
average of some observable A in the state \ip) by (A)^. For At much smaller 
than experimental time resolution, standard quantum mechanics average is 
already an average over e values. Indeed, repeating the measurement on the 
same state, which dehnes a quantum measurement, we also sample e proba¬ 
bility distribution. Actually, in quantum mechanics averaging over repeated 
experiments is the same as an average over simultaneous experiments on a 
large number of identical copies of the apparatus. In our case the two aver¬ 
ages are different in principle, although they coincide for practical purposes, 
unless experiments are performed in time intervals shorter than At. We con¬ 
tinue nevertheless, to use the double averaging notation, to underline the 
conceptual différence. 

Suppose we préparé the System at initial time t = 0. lî h fluctuâtes, po¬ 
sition measurements reçoive additional contributions when averaging over e 
distribution. For mean position value this term is ~ 1 -|- O (a), which 

is time independent, since fluctuation averages are time translation invariant. 
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The same results holds at some final time t, again by time translation invari¬ 
ance. This implies that measurements of the mean distance travelled by a 
particle is insensitive to this order a factor, which cancels in the différence 
(x)^(f) — (x)^(O). The observable effect of averaging over £ is only contained 
in the e term in the state time évolution, see équation (13). Same holds for 
position uncertainty and momentum measurements. 

We choose for reference a gaussian profile peaked at p and variance 




^ — \-ipox/a-iPo{t+fo e{t') dt') / {2mh) 

y^(7rh2)i/4 


Using équations ( |3p| , the mean distance travelled by a particle and uncer¬ 
tainty read 


(x)^(t) - (x)^(O) 
(Ax)J(f) - (Ax)J(O) 



_-L - _ —Tf 

2m2 ^ m2 


(17) 

(18) 


Squared uncertainty displays a new contribution, linear in time. This is much 
alike a Brownian motion with diffusion coefficient 


D 


f + 52/2 

2m? 


(19) 


For ô P, D can be regarded as due to scatterings with mean free path 
{p/m)T. 

First term in r.h.s. of équation (18) is typically expected to be dominant 
over the genuine new effect, but considering massless particles, as photons, 
the linear dispersion relation leads to no wave packet spread. A gaus¬ 

sian wave packet reads in this case 




(7r52)l/4 


( 20 ) 


and to leading (non trivial) order in the expansion of \/l + e{t) one finds 


(x)^(t) - (x)^(O) 
(Ax)2 (f) - (Ax)^(O) 

8 


et , 

C^Tt 


( 21 ) 

( 22 ) 
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This random walk behavior leads to interesting effects in interférence ex- 
periments. Consider a light wave packet impinging a plate pierced by two 
slits, and then observed on a screen, producing an interférence pattern. Since 
waves are detected at some fixed distance L from the plate, the h stochas- 
tic nature can be effectively viewed as a change ôt of the travel time, with 


variance, from équation (22) 


4 4c 


(23) 


with t = L/c the time mean value, neglecting the distance between the slits 
compared to L. Assuming a plane wave with frequency u, light intensity I 
say, at the point on the screen équidistant from the two slits is 


1 

4 


/ oc - |e-*‘^h+<5o) ^ g 


-*a.(t+5t2)|2 = 1 (1 + COS [u{ôh - Ôh)]) , (24) 


where the two contributions corne from the two slits and ôti ^2 are their (un- 
correlated) time shift along the path from slits to the screen. In the standard 
case the two waves show a constructive interférence. Here, averaging over 


ôti^ 2 , from équation (23) 


1 / / u^tL 

/oc- l+exp - — 


(25) 


The interférence term decays exponentially for large L, and asymptotically 
intensity behaves as the two waves were not interfering. This is a genuine 
effect of Planck constant fuzziness, which destroys wave cohérence on dis¬ 
tances u}‘^tL/{4:c) > 1. Reasoning in a similar way one finds that the whole 
interférence pattern changes. 


Notice that, using équation (22), the effect can be cast in terms of a 
fluctuating light speed with variance (Ac)^ 


(Ac) 


CT 

4L ’ 


(26) 


Of course, this holds for relativistic particles only, but not in general. 


The resuit of équation (25) could be tested in long-baseline interferomet- 

aimed at detecting gravitational 


rie experiments, like Virgo or Ligo 24,25 


waves from astrophysical sources. To this end, ail noises, mechanical, seismic 
etc., are kept under an exquisite control, sensitivity being eventually limited 
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in the relevant frequency range by the irreducible shot noise contribution, 
due to photon number N (or wave packet phase) Poisson fluctuations in 
light bunches, AN = \/N. For a small r, équation (25) leads to 


AJ 


,T L 


" 5 “ - 7 “ ’ 

I 4 c 


(27) 


which assuming no détection of r related effects, should be smaller than shot 
noise (sn) 

AI\ _ AN 1 IhuAu 

Ws„“ 


nT L 

< 
4 c 


N ^ y I 

with U and Au the light frequency and bunch bandwidth, or 


(28) 


r < 4 ■ 10“^^ ^ 


km /lO^^HzX 


3/2 


U 


-) 


Au 


10 W 


(29) 


We normalized to visible light frequency and status of the art values for Au, 
L and power I. This bound would translate into a lower limit for an effective 
energy scale of h dynamics A = h/r > 10^° GeV. 


3.2 Harmonie oscillator and cohérent light in cavities 

For a one-dimensional harmonie oscillator with frequency u and mass m the 
Hamiltonian reads 


XX 1 \ 2 , 2 

H = ^(1 + e)Po + —(1 + £)x„. 


(30) 


i.e. a standard harmonie oscillator with time dépendent mass and frequency, 


M = m/(l + e) and G = a;(l + e), IM, 27 , with Mil = mu a constant 


The Hamiltonian dépends on time via an overall multiplicative factor, and 
the commutator [H(t),H(F)] vanishes, implying that Dyson sériés for time 
évolution operator can be computed explicitly 


U{t) = exp (—- J H{t')dt' 


Deflning, with standard normalization, création operator 

1 


(31) 


a = 


mu 


Xo + ï 


\/2mhu 


Po 


(32) 
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we hâve 


H = hu:{l + e) ( a'^a + - 


(33) 


Notice that a and a^ do not explicitly dépend on time and obey 

[a, a^] = 1. (34) 

The a eqnation of motion 

da(t) 


dt 


= -iüo (1 + e{t)) a(t), 


lias the formai solntion 

a(t) = a(0)e“*‘^* 


E 


-lOJ) 


n\ 


dti.... / dtne{ti)....e{tn) ■ 


(35) 


(36) 


Averaging over e{t) probability distribntion and compnting n-point corréla¬ 
tion fonctions in terms of two-point corrélation à la Wick 


a(i) = a(0)e-“‘5^ 




2\k 


2k\ 


i2k-l)\\{TtY 


E^s^ = »(•))' 


—i ujt^—uj^Ttf2 


(37) 


Apart from standard oscillatory term, evolntion is exponentially damped on 
time-scales larger than the characteristic time 2(a;^r)“^. 

Consider now a cohérent state |A) at time t = 0. With no loss of generality 
we take A real. As discussed already, position/momentum measnrements at 


some particnlar time amonnts to measnre xq and po- From eqnation (36) 
after averaging over e: distribntion 


(x)aW = 


2h 


mu 


A cos(a;t)e 


—üj‘^Ttl2 


(p);^(t) = —y/2hmu X sm{ut)e ^ 


(x2) (t) = — 

mu 


(p^)aW = 


^ + A^ (^1 + cos(2u;t)e-"'"*) 
^ + A^ (l - cos(2cnt)e-"^'"*) 


(38) 
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so that cohérent states do not saturate the lower lirait of Heisenberg relation 
as time flows, unless A = 0, 


-r-— h 

AxaAp;, = 2 


1 + 2A2 (l - 


(39) 


This effect is in ail similar to decoherence processes affecting the oscillator 
phase but leaving its energy unperturbed. We show position mean value, 
squared mean value and uncertainty in Figure 1. The growing behaviour of 
Ax^a can be appreciated provided the state remains in a cohérent conûgu- 
ration for times tu > {ru)~^. 

These features could be constrained using optical cavities. Excited by 
an external cohérent beam, they can store cohérence properties of electro- 
magnetic held for long times, and consist in spatial conûnement between two 
highly reflecting surfaces of a well dehned propagation mode. Commonly 
used in laser/maser physics and in optical experiments where high spectral 
spatial purities are required (high resolution spectroscopy, interferometry, 
quantum optics, etc.), their conûnement ability can be quantified, as for me- 
chanical oscillators, in terms of quality factor Q = utc, the ratio between 
energy lost in a cycle to the energy stored in the cavity, with u the frequency 
and te the cavity decay time. Using supermirrors, at optical frequency, values 
of Q ~ 10^® are accessible 


A single mode cohérent state of an electromagnetic harmonie oscillator 
is a minimum uncertainty states for any pair of orthogonal field quadrature 
operators 29,30], the analogue of position and momentum operators. In 


phase space this state is represented by a two dimensional Gaussian dis¬ 
tribution with equal variances at ail direction. To keep a stricter analogy 
with the harmonie oscillator discussed here, we define field quadratures as 
X= -|- al) and Y = i^Jh/2 (al — a). Measurements of the uncer¬ 

tainty for a given quadrature of an electromagnetic mode can be obtained by 
a homodyne detector. 

For a random h, cohérent configuration of radiation do not saturate the 
lower bound h/2 for AX AY, which monotonically increases and is related to 
|Ap, the mean photon occupation number, see équation (39). Measurements 
in résonant cavities however, are limited by te- For t > t^ the cohérent 
electromagnetic field escapes from cavity due to unavoidable couplings to 
the external thermal bath, and the System evolves towards vacuum state 
A = 0. To account for this effect we modify équation (39) by introducing an 
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CO t 


Figure 1: The time évolution of position mean value (x)^ (solid), (x^)^ 
(long-dashed) and squared uncertainty Ax.^x(short-dashed) for a cohérent 
State with A = 1. Values are in units of appropriate powers of the length 
unit x/h/Jjnüj). We hâve chosen an unrealistic large value ur = 0.05 to 
emphasize the non standard time behavior with respect to the case of a 
constant h. 
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exponentially damped A 


AXxAY. = - 


2 L 


l + 2A2e-2'/*= 


(40) 


This approximation is satisfied ii te = Q/OJ is much larger than r, so that A 
adiabatically decays on r time scales. The r.h.s. of équation (40) grows for 
small times, reaches a maximum at 


axaAyJ, =5 


1 + X^Qut 


2 \ {2+Qüjt) / (Qüjt) 

2 + Qut ) 


(41) 


and eventually decays towards the standard value h/2. For Qut <C 1, t* — 
tc/2 - Q'^ t/8. 

Consider one measures the uncertainty product with some error. In the 
standard approach, this product for a cohérent state gives h/2, so is a way 
to détermine Planck constant, with some uncertainty, h ± Ah. If the time 
behaviour of équation (40) is undetected, in particular the peak at t*, this 
means that r should be sufhciently small 

2 


\ Qut 


{2+Qu)t) / (Qut) 

< 


(42) 


2 + Qut J h ’ 

To hâve an order of magnitude of this bound, we take A = 1. Choosing 
Q = 10^®, a; ~ 3 ■ 10^® Hz, 28 , and a measurement uncertainty Ah/h ~ 1% 
we obtain 

r<10"^2s, (43) 

or, in terms of energy scale, A = h/r > 10® GeV. For smaller Ah, the bounds 
on T scales approximately as {Ah/h)/{Qu), see équation (42). 


4 Concluding remarks 

We hâve considered the possibility that Planck constant is a randomly fluc- 
tuating quantity. These fluctuations are viewed as the manifestation at low 
energies of some fuzzy structure of spacetime and phase space at very high 
energy scales E > A, and are found to change the standard results of quan¬ 
tum mechanics, such as the Heisenberg uncertainty relation and interférence 
of massive or massless particle beams. The novel effects can be effectively 


14 
















regarded as due to the coupling of a particular System, such as an harmonie 
oscillator or a free partiale, with some dynamical field. Even if this field has 
reached a stable background configuration, we expect its value to undergo 
small quantum fluctuations, which introduce a noise in the time évolution of 
the System. In particular, we hâve studied the case in which this noise shows 
up in terms of a random Planck constant with a variance times typical cor¬ 
rélation time parameter r = Interestingly, the model can be strongly 

constrained by future experiments with long baseline interferometers or on 
cohérent light cavities, and adopting présent status of the art parameters 
for these kind of experiments, we found that r could be bound in dedicated 
measurements to a very tiny value, of order of 10“^^ — 10“^^ s. These time 
values translate into energy scales A ~ h/r larger than 10® —10^° GeV, which 
are well below the Planck mass scale, 10^® GeV, yet above the energy range 
which can be explored with présent accelerators, like LHG at GERN, and 
presnmably also unaccessible to next génération accelerator programs. 

Of conrse, since Planck constant is nbiqnitous, there is potentially a very 
large nnmber of different dynamical Systems and experimental approaches 
which may be exploited, and which were not covered here. Possibly, some of 
them may be even more powerful than those we discussed in constraining a 
stochastic h. To stndy them ail would entail the reconsideration of the basics 
of ail qnantnm phenomena, an enterprise which is well beyond the scope of 
this paper. 
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